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We study the phase diagram of the (7(2) x (7(2) scalar model in d = 4 dimensions. We find that the phase 
transition is of first order in most of the parameter space. The theory can still be relevant to continuum physics 
(as an effective theory) provided the transition is sufficiently weakly first order. This places restrictions on the 
allowed coupling constants. 



1. The U{N) X U{N) model 

We consider a scalar field theory described by 
the action 
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in Euclidean space, where (t>{x) is a complex 
N X N matrix. The action (|^) is invariant under 
a U{N)l X U{N)ii (to be taken N ^ 2) sym- 
metry, under the global symmetry transforma- 
tion (p L(f)R\ where L,R are U{N) matrices. 
This model has been considered as a low-energy 
effective theory to describe the strong-coupling 
extended technicolor models and top-condensate 
models of electroweak symmetry breaking . 

This model, in contrast to the 0{N) model, is 
known to possess for A2 ^ a first order phase 
transition whose strength varies in the (Ai,A2) 
space. This is due to the well-known Coleman- 
Weinberg instability |^ , and it places restrictions 
on the allowed parameter space of couplings. As 
one adjusts m? past a critical value, m^, the vac- 
uum expectation value (v.e.v.) v jumps discontin- 
uously from zero in the unbroken phase to some 
finite nonzero value in the broken phase. Then, 
if the model is to be a valid low-energy effective 
theory, relevant to continuum physics, the cou- 
plings (Ai,A2), should belong to a region where 
the phase transition is sufficiently weak first or- 
der. It is only then that v can be small compared 
to the cut-off A. 

The model (0) has been studied in perturba- 
tion theory, in terms of the effective potential and 
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the renormalization group (RG) in ref. [| and 
in the language of RG flows and its associated 
fixed points |^J^ in A preliminary investi- 

gation on the lattice was undertaken in ref. 0. 
We present here a more complete investigation of 
the phase structure of the C7(2) x U{2) model, by 
performing Monte Carlo simulations of the lattice 
regularized version of the action (^ above. 

2. Perturbation theory 

In the sequel we take = 2, so the action (|^) 
depends on eight degrees of freedom. More de- 
tails about the model can be found in [^]. If 
A2 = 0, then the symmetry is enhanced to 0(8) 
(0(27V^) for general N). The pattern of symme- 
try breaking depends on the sign of A2 . If A2 > 
then the breaking occurs according to 



U{2)l X U(2)r U(2)v 
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resulting in four Goldstone bosons, while if A2 < 
0, the symmetry breaking pattern is that of 



U{2)l X U(2)r U(lf 
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resulting in five Goldstone bosons. 

Fig. (|l|) displays the RG fiows, at one- loop 
(solid lines) and two-loop (dash lines) level, 
within renormalized perturbation theory. Then, 
starting from bare couplings (Ai(A), A2(A)), if 
A2 ^ all RG fiows in the infrared intersect 
the "stability line" eventually becoming run- 
away trajectories. The phase transition is then of 
first order ||^. For A2 small, though, the fiow is 
rather slow and even though the "stability line" is 
crossed, this happens after many decades of run- 
ning; the phase transition is, in this case, weakly 
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Figure 1. Perturbative RG trajectories starting 
from bare couplings (Ai(A), A2(A)) along the lines 
Ai = 2 or A2 = 2. The solid lines (dash lines) cor- 
respond to one-loop (two-loop) trajectories while 
the stability line is indicated as a dotted line. In- 
dicatively, the dots along a trajectory represent 
the evolution of couplings after running by a fac- 
tor of e down to the infrared. 

first order, with v ^ A. Along the A2 = axis 
the transition is known to be second order, and 
the well known triviality analysis ^ applies. The 
two-loop corrections seem to improve the hierar- 
chy A/v, as found also in In particular, there 
exists a region with Ai, A2 > where the flow is 
towards larger values and it appears that it never 
crosses the stability line. However, this only hints 
upon the breakdown of perturbation theory and 
a nonperturbative analysis is called for. It is also 
essential that no other fixed point, unreachable 
in perturbation theory, exists. Should one be 
present, the RG trajectories would be distorted 
and there could be regions where the transition 
is second order. We found no evidence of such a 
fixed point. 

On the lattice, the physical parameter control- 
ling the running of the couplings and hence the 
size of the corrections is the correlation length ^ 
of the system. Then one expects that A/v ~ £^^, 
where /? is the appropriate critical exponent: a 
large hierarchy will only be possible if ^ is big, 
or equivalently, that the transition is weakly first 
order or second order. 



Table 1 

Estimates of the jump in the order parameter 
evaluated at the critical point m^, and of the cor- 
relation length estimated by ^ ~ L* /2, where 
L* is the smallest lattice where coexistence was 
found, or from the effective potential. 
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(0.5, -0.45) 


-0.772 


0.83 
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(-0.22, 0.5) 


-0.91 


2.10 
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(0, 0.5) 


-2.42 


0.5 


40 


(-3.97, 8) 


-1.30 


10-20 


< 2 


(-14.97, 30) 


-1.50 


20-40 


< 2 


(0, 8) 


-24.75 


0.15 
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(0, 16) 


-43.98 


0.15 
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(0, 30) 


-77.0 


0.19 
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(8, 8) 


-63.8 


0.11 
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(8, 16) 


-82.5 


0.16 


6 


(8, 30) 


-114.5 


0.17 
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3. Monte Carlo Results 

Table I shows all the points in the (bare) 
coupling constant space at which we performed 
Monte Carlo simulations. We used two differ- 
ent programs checked against each other: one 
based on a simple one-hit Metropolis algorithm 
with a uniform step tuned so that the acceptance 
rate is about 60% and the other based on the 
hybrid algorithm (with or without Fourier accel- 
eration). In this latter method, two parameters 
have to be chosen, namely the number of leap- 
frog steps and the step size. This allows better 
control of the autocorrelation time. We found 
optimal CPU performance for 5-8 leap-frog steps 
before each Metropolis test. The hybrid one per- 
formed clearly better. 

We used as an order parameter (following 0|) 
the expectation value of the U{N) x U{N) invari- 
ant operator O = Tr^^^ which corresponds to 
the susceptibility, where (f> is the lattice average 
of each field component. < O > is proportional to 
in the broken phase and zero in the unbroken 
phase, modulo finite size corrections. 

We used lattices of sizes ranging from = 4'' 
to = 14^. In order to obtain information about 
the order of the transition, at each given (Ai, A2) 
we searched for hysteresis effects in the measure- 
ment of the order parameter by performing ther- 
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mal cycles in the relevant parameter, m^, across 
the critical region. Strong hysteresis loops is an 
indication of a strong first order transition. 

On smaller lattices, (4*, 6*, 8''), the critical re- 
gion was identified by searching for a double-peak 
signal in the histogram distribution of Tr(0(/)^). 
We then moved to bigger lattices (10^, 12'*, 14*) 
to look for coexistence. Along the process of in- 
creasing the lattice size we eventually begin to see 
metastability at some size L* . We estimate then 
the correlation length to be ^ ^ L* /2. Crude as 
this procedure may seem, it is physically mean- 
ingful and it agrees, where comparison is possible, 
with the effective potential. 

Our results are as follows. All points close 
to the stability line exhibited marked hysteresis 
loops and hence show strong first order transi- 
tions, becoming stronger as we move up along 
the stability line. For such couplings, ^ ~ 1 so 
the cut-off effects are big and the connection to 
continuum physics questionable. In the weak cou- 
pling region (Ai,A2 < 1), we also computed the 
one-loop bare effective potential Q and found 
that it agreed with the numerical data within 
10 — 30%. Next we investigated couplings along 
the A2 axis. Typically, runs on 4'',6* lattices 
did not show any hysteresis effects. However we 
found clear sign of the existence of two minima on 
lO'* — 14''' lattices. Fig. (Q) displays tunneling be- 
tween two phases for the point (0, 8) on a 10* lat- 
tice. For the same point, clear signal was found of 
two coexistent minima on the 12* lattice, but no 
tunneling was observed. The transition becomes 
stronger with increasing A2. For points deep in 
the Ai > 0, A2 > region, we were able to ob- 
serve coexistence of phases, but only on 12*, 14* 
lattices. The transition is always clearly first or- 
der, but characterized by correlation lengths, as 
expected, larger than those obtained close to the 
stability line. For couplings close to the Ai axis, 
our numerical results are in agreement with the 
expectation of the RG for a weak first order tran- 
sition. In this case, though, it is difficult to dis- 
tinguish a weak first order from a second order 
transition. 

Our results are summarized in Table ([T]) . From 
these, one can see that in most of parameter space 
the v.e.v., u(A), is typically only one order of 
magnitude smaller than the cut-off. Our results 
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Figure 2. Order parameter time history (over 2 x 
10^ hybrid sweeps) for the point (0, 8) on a 10* 
lattice at = -24.69. 



are consistent with the standard perturbative pic- 
ture of first order phase transitions and the ab- 
sence of any nontrivial fixed point. The hierar- 
chy A/w is not "tunable" by m? as in the 0{N) 
model, but rather depends on A2 . Phenomenolog- 
ically viable models must lead to couplings with 
small A2 in order to support a large hierarchy. 
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